caixeta <- read.csv("caixeta.csv", header=TRUE, as.is=TRUE)
head(caixeta)
tail(caixeta)
summary(caixeta)
dim(caixeta)
class(caixeta)
names(caixeta)

library(graphics)

head(caixeta$cap)

caixeta$dap <- (pi/4)* (caixeta$cap/10)
head(caixeta$dap)
length(caixeta$dap)

caixeta.2 <- aggregate(caixeta$dap,by=list(local=caixeta$local,parcela=caixeta$parcela,arvore=caixeta$arvore,especie=caixeta$especie),FUN=sum)

dim(caixeta.2)
head(caixeta.2)

names(caixeta.2)[5] <- "dap.total"
caixeta.2

str(caixeta.2)

h.media <- aggregate(caixeta$h,by=list(local=caixeta$local, parcela=caixeta$parcela,arvore=caixeta$arvore,especie=caixeta$especie),FUN=mean)
head(h.media)

sum(caixeta.2[,1:4]!=h.media[,1:4])

caixeta.2$h <- h.media$x

head(caixeta.2)

unique(caixeta.2$local)

caixeta.tabebuia <- subset(caixeta.2, especie=="Tabebuia cassinoides", select =c(local, parcela, arvore, especie, dap.total, h))
caixeta.tabebuia

plot(dap.total~h, data=caixeta.tabebuia, subset=local=="retiro")
scatter.smooth(caixeta.tabebuia$dap.total~caixeta.tabebuia$h, subset=caixeta.tabebuia$local=="retiro")

plot(dap.total~h, data=caixeta.tabebuia, subset=local=="jureia")
scatter.smooth(caixeta.tabebuia$dap.total~caixeta.tabebuia$h, subset=caixeta.tabebuia$local=="jureia")

plot(dap.total~h, data=caixeta.tabebuia, subset=local=="chauas")
scatter.smooth(caixeta.tabebuia$dap.total~caixeta.tabebuia$h, subset=caixeta.tabebuia$local=="chauas")

library(lattice)

xyplot(dap.total~h|local, data=caixeta.tabebuia)
